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Helium nanodroplets and trapped Bose-Einstein condensates in dilute atomic gases offer com- 
plementary views of fundamental aspects of quantum many-body systems. We discuss analogies 
and differences, stressing their common theoretical background and peculiar features. We briefly 
review some relevant concepts, such as the meaning of superfluidity in flnite systems, the behavior of 
elementary excitations and collective modes, as well as rotational properties and quantized vorticity. 



I. INTRODUCTION 



^ I Helium clusters have been the object of a rather extensive investigation in the last two decades. They are becoming 

even more interesting nowadays, due to the improved accuracy of recent experiments where several properties of 
CN ' the clusters are characterized by using dopant atoms and molecules as efficient probes M. In 1995, the observation 
of Bose-Einstein condensation (BEC) in ultracold vapors of alkah-metal atoms |Q-Q opened another active field of 
research. This event is seemingly disconnected from the physics of helium clusters. The framework in which BEC 
was initially investigated is a mix of atomic physics and quantum optics, while helium clusters are a typical subject 
fvi ' of condensed matter and chemical physics. As a consequence, some of the key words used to interpret the first 
f^ ] experimental results on BEC (like order parameter, coherence length, phase coherence, elementary excitations, and 
f~^ ' others) had sometimes different meaning in the two communities. However, as soon as trapped condensates were 
^^ ■ sufficiently characterized, by using the Gross-Pitaevskii theory as a starting point, most of the physics behind these 
systems appeared to be strongly connected with the physics of superfluids. One can even say that helium droplets 
and trapped atomic gases correspond to two limiting cases of the same system since they are, respectively, examples 
of very dense and dilute finite-sized quantum fluids. Despite the different type of confinement (self-binding for helium 
droplets and external trapping for atomic gases) these systems belong to the same conceptual framework, possibly 
based on a common language. 

i-rt ■ A quantum fluid is a system of many interacting particles for which one can locally define both a density and a 

r^ ' velocity field, whose behavior is strongly affected by quantum correlations. An example of these quantum effects is the 

O ■ zero point motion, which is crucial in determining the typical lengthscales for density modulations, like the surface 

J-^ ' thickness at the fluid boundaries, the core size of a vortex, the size of soliton-like structures. Another important 

effect is the fact that the fluid, or part of it, moves with a non-dissipative and irrotational velocity field, which means 

that the system is superfluid. These and other features can be associated with the existence of an order parameter. 

Qualitatively, as the particles are cooled down to enough low temperatures their de Broglie wavelength becomes larger 

j^ than their average separation, so that they loose their individual identity. In the case of bosons, they can be treated 

as components of a single "macroscopic wave function" . The concept of order parameter is the appropriate tool to 

implement this qualitative idea. 

The temperature at which these effects occur is generally very low, so that the fluid actually freezes before reaching 
this regime. In the case of liquid ^He, the light atomic mass allows the system to remain liquid down to zero 
temperature and the system becomes a quantum fluid for temperatures lower than about 2.2 K. In the case of trapped 
condensates, vapors of alkali-metal atoms are kept in a metastable gas-phase down to temperatures of the order of tens 
of nanokelvin, where quantum effects becomes dominant. The fact that these systems are in a metastable conflguration 
is compatible with the existence of a kinetic equilibrium which is ensured by two-body elastic collisions. The transition 
to the thermodynamic equilibrium, given by the crystal phase, is instead driven by three-body recombinations, which 
are however rare events in very dilute samples and take place on a relatively long time scale. 

When the fluid becomes a quantum fluid the type of quantum statistics is also crucial and one can flnd either 
a Fermi or a Bose fluid with very different properties. Liquid "^He is the fermionic counterpart of *He, and also 
alkali-metal atoms with both Fermi and Bose statistics can be trapped and cooled down to the quantum regime. In 
this brief review, however, we concentrate mainly on Bose fluids. 

The paper is organized as follows. First we remind the basic notions of Bose-Einstein condensation, as a common 
starting point. Then we discuss two opposite limits: a dilute and cold gas, where one can develop a rigorous theory for 
interacting bosons, and a dense quantum fluid, where more sophisticated many-body techniques are required. Once 



the theoretical framework is sketched, we finally describe some predicted and/or observed properties of both helium 
droplets and trapped condensates (critical temperature, ground state, excitations, superfluid effects, etc.) in order to 
emphasize relevant analogies and differences. 

II. BOSE-EINSTEIN CONDENSATION 

The intimate link between helium droplets and trapped condensates is the occurrence of Bose-Einstein condensation 
and the fact that this phenomenon plays a key role in determining most of the properties of both systems. 

In the case of noninteracting particles, Bose-Einstein condensation has well established features. The system can 
be represented by a set of single-particle states of energy e^, each one occupied by a certain number of particles rii. 
At a given temperature T, the average occupation numbers are fixed by Bose statistics, resulting in the well known 
law (ni) = l/{exp[/3(ei — /i) — 1]}, where (3 = (kBT)~^ and /i is the chemical potential. The latter is fixed by the 
conditions that the total number of particles, N ^ is equal to the sum J^ti^i)- F^^ relatively high T one finds that 
this condition is satisfied for /i < eo, where eo is the energy of the lowest single-particle state, and the particles are 
thermally distributed over many states, all of them having (n^) of the order of 1 (or, equivalently, much less than N). 
However, there is no exclusion principle for bosons: many of them can stay in the same state. This happens when 
the temperature is lowered; when the chemical potential approaches eo, the occupation of the lowest state becomes a 
macroscopic number Nq = (uq) of the order of N. This phenomenon is called Bose-Einstein condensation. The critical 
temperature at which it occurs can be calculated starting from the knowledge of the single-particle Hamiltonian and 
keeping the appropriate thermodynamic limit TV — > oo. 

In the textbook case of a uniform system of free bosons of mass m, one increases the volume V while keeping the 
density n — N/V constant; in this way, one gets the critical temperature g 

''^' ^ ^ Kwm) ^'^ 

where C, is the Riemann zeta function and C(3/2) ~ 2.61. This result can be rewritten as nAi^ ~ 2.61, where 
At = [27r?i^/(mfcBT)]^/^ is the thermal de Broglie wavelength. This is a quantitative translation of the qualitative 
idea mentioned in the introduction: quantum effects become crucial when the de Broglie wavelength becomes larger 
than the average separation between particles. The same calculations also yields the condensate fraction (Nq/N) = 
1 — [T /Tcf'/'^ . At T = all particles are in the eo = state, that is, the system is fully condensed. This state has a 
uniform density and is characterized by having zero momentum; for this reason, one often speaks of condensation in 
momentum space. It is important to stress here that, when applied to liquid ^He {m and n being the atomic mass and 
the liquid density, respectively), Eq. (nf) gives a transition temperature reasonably close to T\, i.e., the temperature 
at which helium becomes superfluid. This very simple and striking result marked the beginning of the theory of 
superfluids, starting from the pioneering work of Fritz London [pj. 

In the case of N bosons confined in a spherical harmonic potential V^xt = (l/2)'Tia;^Qr^ , the proper thermodynamic 
limit is obtained by letting N ^r oo and Who — > 0, while keeping the product Nuj"^^ constant. In this way the critical 
temperature is well defined; one finds 

keT, = n^ho [tF^) - 0-94 ?Jc^ho N^'^ (2) 

and the condensate fraction is {Nq/N) = 1 — (T/Tc)^. Notice that the different T-dependence exhibited by the 
condensate fraction compared to the uniform case is the consequence of the higher density of states characterizing the 
harmonic oscillator Hamiltonian. At T = all particles are in the lowest eigenstate of the harmonic oscillator, namely, 
/i ^ eo = (3/2)?i(jJho- Thus the particle density has the form of a Gaussian and the same is true for the momentum 
distribution. Differently from the uniform gas, condensation occurs in both real and momentum space. Actually, the 
occurrence of a sharp peak in the velocity distribution of a cloud of '^^Rb atoms, released from a magnetic trap, was 
the first evidence of Bose-Einstein condensations in these cold gases M , and the transition temperature was found to 
be very close to the value predicted by Eq. (Q). 

Since both liquid "^He and dilute atomic gases are systems of interacting particles, a crucial question concerns the 
role played by the interaction, that is, whether and how much the interatomic forces modify the properties of Bose- 
Einstein condensation. Understanding the interplay between quantum statistical and dynamic effects represents the 
challenging task to achieve both experimentally and theoretically in these quantum many-body systems. 



Let us consider a system of N particles interacting via a potential ^(r — r') and, possibly, with an external field 
T4xt(r). In principle, one should solve the exact many-body Schrodinger equation or, equivalently, find the eigenvalues 
of the Hamiltonian 

/^2 "I 1 /" 

dr*^(r) V^ + ycxtfr) "^(r) + - dr'dr ¥{r)iiUr')V(r - r'W(r')i/(r) (3) 

2m J 2 J 

where ^(r) and ^'l'(j-) are the boson field operators that annihilate and create a particle at the position r, respectively. 
This is a difficult task and we do not discuss here the possible strategies. We only want to emphasize the meaning of 
Bose-Einstein condensation in terms of a generic solution of the many-body problem. An elegant way g consists in 
writing the one-body density matrix, which is defined as 

p(r',r) = (*t(r')^(r)) (4) 

where the average is taken in the state that one wants to describe. This quantity characterizes the correlations between 
particles located in different points in space. One can use a complete set of single particle states to explicitly write 
p(j' , r) as a matrix. Then one can diagonalize it, finding eigenstates and eigenvalues. The latter can be interpreted as 
the occupation numbers associated with each eigenstate. In analogy with the ideal Bose gas, one might expect that 
below some critical temperature one of these eigenvalues, instead of being of order 1 as all the others, is a number 
Nq of the order of N. The corresponding eigenstate then has a special role in the system, since it is macroscopically 
occupied. When an interacting system exhibits such a behavior, one says that it is Bose-Einstein condensed. 

In the case of a uniform system, the appropriate basis is a set of plane waves, eigenstates of the momentum p. If 
Nq is the number of particles condensed into the lowest energy p = state, then the density matrix can be written 
in the form 

p(r',r)=7Vo+p(r-r') (5) 

with /5(r — r') — ^ Up exp[ip • (r — r')/h]. This summation is made over all states except p = 0; for r — r' = it 
simply gives the total number of particles with p 7^ 0, that is p{0) — N — Nq. In the opposite limit, |r — r'| -^ 00, due 
to destructive interference between the various phase factors, the function j5 vanishes, while the condensed part of the 
density matrix remains everywhere constant. This infinite range of the the density matrix (off-diagonal long range 
order) is also called first order coherence. It is a key consequence of Bose statistics and is the origin of the quantum 
correlations that are responsible for superfiuidity. This is the basic physics behind the superfluid behavior of liquid 
''He below T\. In such a dense system, the effect of the interaction, which is implicitly included in Eq. (^, is very 
important. Compared to the ideal gas, it changes the features of the phase transition and gives a large depletion of 
the condensate, Nq/N being of the order of 10% even at T = 0. 

This description can be easily generalized to nonuniform systems. In this case, the contribution to p{r' , r) associated 
with the condensed state is no more a number, but depends on positions. The straightforward generalization of Eq. (|5|) 
is 

p(r',r) ==$*(r)$(r')+p(r,r') (6) 

where <I>*(r)<I>(r') is of order N, while p{r,r') vanishes for large |r— r'| g]. This is equivalent to say that the bosonic 
field operator splits in two parts 

*(r) = $(r) + *'(r) , (7) 

where the first term is a complex function associated with the macroscopic occupation of the condensate and the 
second one is the field operator associated with the noncondensed particles. The complex function can be defined as 
the expectation value of the field operator: 

(#) = $ = |$|e'^ . (8) 

It behaves as a classical field having the meaning of an order parameter, and is often named macroscopic wave function. 
Its modulus gives the condensate density through rio(r) — |<I>(r)p, while the phase S can be used to define a velocity 
field through v = {h/m)W S. 

Similarly to the case of uniform gases, the fact that the order parameter has a well-defined phase corresponds to 
assuming the occurrence of a broken gauge symmetry in the many-body system. It is worth noticing that, strictly 
speaking, in a finite-sized system neither the concept of broken gauge symmetry, nor the one of off-diagonal long-range 
order can be applied. In fact, the density matrix (g) vanishes for distances of the order of the size of the system so 



that, in order to make the separation of the condensate and noncondensate components physically meaningful, it is 
crucial that the lengthscale at which p(r,r') vanishes be much smaller than the size of the system. The condensate 
wave function $ can nevertheless be always calculated as the eigenfunction of the density matrix with the largest 
eigenvalue. A numerical diagonalization of the density matrix was performed, for example, by Lewart et al. flOfl in 
order to estimate the condensate density, no(r), in helium droplets with up to 240 atoms, using a variational Monte 
Carlo approach. 



III. DILUTE GASES AND GROSS-PITAEVSKII THEORY 



The decomposition of the field operator in Eq. (R) is particularly useful when the depletion of the condensate, i.e., 
the fraction of noncondensed particles, is very small. This happens when the interaction is weak, but also for particles 
with arbitrary interaction, provided the gas is dilute. In this case, one can expand the Hamiltonian \^ by treating 
^' as a small quantity. 

In a uniform gas the condensate wave function, $, is just a constant and expanding H to the lowest order in *&' 
gives access to the equations for the excited states. This procedure was introduced by Bogoliubov ||lj], who found an 
elegant way to diagonalize the Hamiltonian by using simple linear combinations of particle creation and annihilation 
operators. These are known as Bogoliubov's transformations and stay at the basis of the concept of quasiparticle, 
which is one of the most important concepts in quantum many-body theory. 

If the gas is nonuniform, the theory has to include an equation for the spatial variation of the condensate, even at 
the zeroth-order in ^P'. A possible strategy consists in writing the Heisenberg equation for the evolution of the field 
operators, using the Hamiltonian (j^): 



in^^i>{v,t) = [i>,H] = [- 



2m 



+ Kxt(r)+ / dr' ¥{r',t)V{r'-r)^{r',t) *(r,i) 



(9) 



Then the lowest order is obtained by replacing the operator ^ with the classical field $. In the integral containing the 
atom-atom interaction V{r' — r), this replacement is, in general, a poor approximation when short distances (r' — r) 
are involved. In a dilute and cold gas, one can nevertheless obtain a proper expression for the interaction term by 
observing that, in this case, only binary collisions at low energy are relevant and these collisions are characterized by 
a single parameter, the s-wave scattering length, independently of the details of the two-body potential. This allows 
one to replace V{r' — r) in (Q) with an effective interaction V{r' — r) = gS{r' — r) where the coupling constant g is 
related to the scattering length a through g = Anh a/m. Using this pseudo-potential in Eq. (||) and replacing ^ with 
$, one gets the following closed equation for the order parameter: 



* 9 ^ , 
in—^{v,t) 



2m 



+ Kxt(r)+g|$(r,t)|2 U(r,t) 



(10) 



This is known as Gross-Pitaevskii (GP) equation [|2|,^. It has been derived assuming that N is large while the 
fraction of noncondensed atoms is small. On the one hand, this means that quantum fluctuations of the field operator 
have to be small, which is true when n|ap <C 1, where n is the particle density. In fact, one can show that, at T = 
the depletion of the condensate is proportional to (njap)"'^/^ pj. On the other hand, thermal fluctuations have also 
to be negligible and this means that the theory is limited to temperatures much lower than Tc- Within these limits, 
one can identify the total density n with the condensate density uq. 

The stationary solution of Eq. (jl^) corresponds to the condensate wave function in the ground state. One can write 
<I'(r,t) — (j){r)exp{—i^t/h), where /i is the chemical potential and </> is real and normalized to the total number of 
particles, J dr cj)'^ = Nq = N. Then the GP equation becomes 



2m 



+ V;,t(r)-l-g0'(r)U(r)=Mr) 



(11) 



This has the form of a "nonlinear Schrodinger equation" , the nonlinearity coming from the mean-field term, propor- 
tional to the particle density n(r) — 0^(r). It is worth noticing that the same equation can be obtained by minimizing 
the energy of the system written as a functional of the density: 



E[n] ^ dr 



2m 



|V^AI|2 + nKxt(r) + 



gn^ 



E, 



qp 



E,. 



E\r 



(12) 



The first term corresponds to the quantum kinetic energy coming from the uncertainty principle; it is usually named 
"quantum pressure" and vanishes for uniform systems. 

In the absence of interactions (g = 0), the GP equation (|l^) reduces to the usual single-particle Schrodinger 
equation. Conversely, when g ^ 0, the mean-field contribution can have an important role in determining the energy 
and density distribution of the condensate, despite the smallness of the gas parameter n|ap. In fact, what matters is 
the relative weight of the mean-field potential, the kinetic energy, and the external potential. As we will see later on, 
one can easily find situations where the gas is dilute but strongly nonideal. It is also worth pointing out the key role 
played by the chemical potential in the GP theory, where the time dependence of the stationary order parameter is 
fixed by fj, and not by the energy. This is a consequence of the fact that the order parameter is not a wave function 
and that the GP equation is not a Schrodinger equation in the usual sense of quantum mechanics. From the point of 
view of many-body theory the order parameter corresponds to the matrix element of the field operator between two 
many-body wave functions containing, respectively, N and N -\- 1 particles. This implies that its time dependence is 
fixed by the factor exp{—i[E{N -I- 1) — E{N)]t} and hence by the chemical potential /i = dE/dN rather than by the 
energy E. 

The excited states of the condensate, at T = 0, can be also calculated starting form the time dependent GP equation 
(O), by looking for small deviations around the ground state in the form 

$(r, t) = e"'^*/'* [0(r) + w(r)e-*'^* + w*(r)e'"*] . (13) 

By keeping terms linear in the complex functions u and w, Eq. (lid) becomes 



huju{v) = [Ho-fi + 2g<j>^{r)]u{r) + gc/)^ {r)v{r) (14) 

-hLJvir) = [Ho-n + 2502(r)]«(r) + g4>^{v)u{T) . (15) 

where i?o — — (?i^/2m)V^ + V^xt(r). These coupled equations allow one to calculate the eigenfrequencies lj and hence 
the energies e — huj of the excitations. In a uniform gas, the amplitudes u and v are plane waves and one recovers 
the Bogoliubov's spectrum ||l^ 



where q is the wavevector of the excitations. For large momenta the spectrum coincides with the free-particle energy 
h'^q'^/2m. At low momenta Eq. (|l6|) instead gives the phonon dispersion lj = cq, where c = [gn/m]^/^ is the sound 
velocity. This transition between the two regimes occurs when the typical wavelength of the excited states becomes 
of the order of the so called healing length, 

e = [Snna]-^/^ , (17) 

which is a rather important lengthscale for superfluidity. When the order parameter is forced to vanish at some 
point (by an impurity, a wall, or something else), the healing length is the typical distance over which it recovers its 
bulk value. In a nonuniform condensate the lowest excitations are not plane waves, but they have still a phonon-like 
character, in the sense that they involve a collective motion of the condensate, and the transition from phonon-like 
to single-particle excitations is still an important feature of these systems. 

An important point which deserves to be stressed is the existence of a close connection between the dynamics of a 
condensate, expressed by the GP equation, and the hydrodynamic equations for an irrotational and viscousless fluid. 
The natural variables for the hydrodynamic description are density and velocity field, which can be related to the 
modulus and the phase of the order parameter as in Eq. (H) . Using these definitions into the GP equation (noh and 
neglecting the contribution of the quantum pressure, V^^/n, one gets the two coupled equations 

^n + V • (vn) = (18) 

™^v + VU+^-)=0 (19) 

where ^(r, t) = gn{r, t) + Vcxt(j)- These are the continuity and Euler equations for an irrotational and viscousless fluid 
in the coUisionless regime, that is, the hydrodynamic equations for a superfluid at T = 0. Neglecting the quantum 
pressure term in the GP equations means that Eqs. dlq) and (n9h are valid for macroscopic motions of the condensate, 
i.e., excitations with wavelength larger than the healing length. 



The equations written in this section stay at the basis of most of the theoretical papers recently written on BEC in 
trapped gases (see |14|] for a recent review) . Analytical and numerical methods have been developed in order to solve 
both the stationary and time-dependent Gross-Pitaevskii equation and the results are often in excellent agreement 
with experiments. On the other hand, the accuracy of the approximations made in deriving the GP equation can 
be directly tested by performing Monte Carlo simulations, which are exact within statistical errors. For trapped 
condensates, this has been done in ||l^ using a Path Integral Monte Carlo method, showing that the stationary GP 
equation is very accurate for T less than about 0.75Tc and providing quantitative estimate for the effects of correlations 
beyond mean-field theory. Such corrections has been also calculated for a uniform hard-sphere Bose gas in |16[| using 
a Green Functions Monte Carlo method. 

It is worth mentioning that the theory of weakly interacting Bose gases born well before the observation of BEC 
in trapped gases, as an attempt to explain the behavior of supcrfluid helium. In the case of helium, however, the 
parameter n\a\^ is of the order of 1, so that the theory looses most of its predictive power. 

IV. DENSE FLUIDS AND THE QUANTUM MANY-BODY PROBLEM 

When the system is neither weakly interacting nor dilute, the mean-field theory based on Eq. (|lO| ) is no more 
applicable. Quantum correlations beyond mean-field are essential; they can give, for instance, a large depletion of 
the condensate even at T = 0, so that perturbative schemes based on the smallness of ^' in Eq. (|^) fail. In order to 
overcome this obstacle, one can follow two alternative and complementary paths: i) keep the theory at a microscopic 
level by solving the many-body problem as accurately as possible; ii) develop phenomenological theories, starting from 
a more macroscopic view of the liquid. 

Among the microscopic approaches one can distinguish different strategies depending on the use or non-use of 
stochastic Monte Carlo procedures. Example of non-stochastic approaches are variational methods and perturbative 
schemes. Most of them relies on diagram expansion and summation techniques typical of quantum many-body 
theories. Jastrow-Feenberg [ |l7[ wave functions are often used in this context as a starting point. They can be 
combined with integral equations as in the optimized hypernetted chain (HNC) theory or with perturbative schemes 
as in the correlated basis function (CBF) theory. The energy minimization within a certain class of variational wave 
functions can be performed with stochastic techniques as in the variational Monte Carlo (VMC) method. This allows 
one to include more sophisticated correlations in the trial wave function as, for instance, in the recent Shadow Wave 
Function theory Monte Carlo algorithms can finally be used to get exact results with statistical uncertainties. This 
is the case of approaches like Green's function MC, Diffusion MC and Path Integral MC, which are powerful, reliable, 
and computationally demanding. It is not our aim to present an overview of all these theories and calculations. 
Appropriate references can be found in the various papers where these techniques are applied to helium droplets (see, 
for instance, [Esj and references therein). 

Conversely, a problem that deserves to be mentioned at this point is whether it is possible to make a bridge between 
dilute gases and dense fluids, like helium, within a single theoretical scheme. The recent Variational Monte Carlo 
calculation by DuBois and Clyde p9[ is an interesting effort in this direction. By using a rather simple variational 
wave function, they evaluated the ground state properties of a trapped Bose gas made of hard-sphere of radius a, 
over a wide range of the parameter na'^. In particular, through the diagonalization of the one-body density matrix, 
they calculated the condensate fraction showing how it decreases from 1 to about 10% in bulk, when na^ goes from 
to the value appropriate for liquid helium, which is about 0.2-0.3 (the range of the interaction is about 2-3 A, while 
the density is approximately 0.022A~'^). They also pointed out interesting features about its spatial dependence, i.e., 
the tendency for dense ffuids to have larger condensate fraction at the surface where the density is smaller. This 
enhancement of the condensate was first found in the Variational MC calculations of Ref. iQ in the case of helium 
droplets. 

Phenomenological approaches can also be applied to study dense systems. For confined fiuids the simplest one is 
the liquid drop model. One considers a droplet of constant density within a radius R and a sharp surface. Then 
one can write the energy as a power expansion in R (or, equivalently, in power on the number of particles N) taking 
bulk energy, surface energy, curvature energy, etc., as input parameters. In the limit of large N, one can consider the 
dynamics of these droplets as due to large wavelength, collective motions. Then, it makes sense to use hydrodynamic 
equations with appropriate boundary conditions in order to get the dispersion of bulk (phonon-like) and surface 
(ripplon-like) excitations pG,El|. Several properties of quantum ffuid droplets are well predicted by this type of 
models ||22| ; what is missing is the accurate description of structures and excitations on a more microscopic scale, i.e., 
at the level of atomic distances. A step further in this direction is the density functional approach. 

In density functional theories at zero temperature, the energy of a system is assumed to be a functional of the 
particle density n{r). A good starting point for Bose fluid is given by 



E[n] = Ec[n] + dr 



2m 



(20) 



which is the natural generahzation of the Gross-Pitaevskii functional (12). Ground state configurations are obtained 



by minimizing this energy with respect to the density. This leads to the Hartree-type equation 

(-^ + Kxt(r) + U[n{r)]^ V^) = mv^^ , (21) 

where U[n] = SEc/Sn{r) acts as a mean field, while the chemical potential fj, is introduced in order to ensure the 
proper normalization of the density to a fixed number of particles. Equation_(pl|) exhibits a formal analogy with the 
Gross-Pitaevskii equation (^l|). The main conceptual difference is that Eq. (|2l| ) is an equation for the density and not 
for the order parameter. While in a dilute gas setting n = |$p is an excellent approximation, in a correlated system 
it is meaningless. 

The correlation energy Ec[n\ incorporates the effects of dynamic correlations induced by the interaction. The fact 
that this energy is a unique functional of n(r) is ensured by the Hohcnberg-Kohn theorem |£^. In the case of dilute 
gases, one simply has Ec[n] — J dr{l/2)gn'^, as in Eq. (^. Conversely, since hquid helium is a strongly correlated 
system, a rigorous derivation oi Ec[n\, starting from first principles, is not available. One then resorts to approximate 
schemes for the correlation energy (see, for instance, Refs. i23,p^ for discussions of this problem from the viewpoint 
of microscopic many-body theory). A simple but useful approach consists of writing a phcnomcnological expression 
for the correlation energy, whose parameters are fixed to reproduce known properties of the bulk liquid. A functional 
of this type was introduced in Ref. |^,|3 to investigate properties of the free surface and droplets of both *He and 
■^He. The correlation energy was written as 



E,[n]^ f dr 



^n^ + ^n^+i + d{Vnf 



(22) 



where 6, c and 7 are phenomenological parameters fixed to reproduce the ground state energy, density and com- 
pressibility of the homogeneous liquid at zero pressure, and d is adjusted to the surface tension of the liquid. The 
first two terms correspond to a local density approximation for the correlation energy, while nonlocal effects are in- 
cluded through the gradient correction. Non-locality effects have been included in DFT in a more realistic way by 
Dupont-Roc et al. [g8| , who generalized Eq. ( |2^ ) to account for the finite range of the atom-atom interaction. Further 
improvements has been added in Ref. |2^. These and other similar functionals have been used for calculations of 
ground state properties of liquid helium in different geometries (see pO| and references therein) . 

Dynamical problems can also be faced with density functionals. On the one hand, one can include in the functional 
terms which explicitly depends on the current density, as done in [E9[ , thus obtaining two coupled equations for both 
density and velocity field, which can be numerically solved. Excited states of helium droplets have been obtained in 
this way |3llj3^' the linearized equations of motions, in this case, correspond to a generalization of the Bogoliubov 



equations^IJP(I5), formally equivalent to the equations of the Random Phase Approximation (RPA) Q with an 
effective interaction of phenomenological nature. On the other hand, one can simplify the problem by using a purely 
hydrodynamic approach, in which the equations of motions for the velocity field are the ones of an irrotational and 
viscousless fluid, but with the density taken from the minimization of the density functional ( po[ ) . This viewpoint was 
chosen, for instance, in Ref. |p3[ to predict the moment of inertia of doped helium droplets. 

V. LIKE AND UNLIKE FEATURES OF HELIUM DROPLETS AND TRAPPED CONDENSATES 

A. Temperature 

In typical experiments the temperature of ^He droplets is about 0.37 K |Q, as the result of a spontaneous evapora- 
tion of the warmest atoms after the formation of droplets in the free jet expansion of helium p5[ . This temperature is 
well below the bulk superfluid transition temperature, T\ = 2.17 K. Path Integral Monte Carlo calculations [Q have 
shown that even small droplets, with few tens of atoms, are superfluid at this temperature. Moreover, the droplet 
is cold enough to neglect the thermal activation of excited states. The lowest ones are discretized collective modes 
(bulk and surface oscillations) which are not significantly populated at T ~ 0.4 K. Thus, for most purposes, helium 
droplets can be considered as T = systems. 

Trapped condensates can be cooled down to few tens of nK. Again this is the result of evaporation, but here the 
evaporative cooling process is induced and controlled through an external rf-field which acts on the shape of the 



confining potential, lowering the edge of the trap and letting the hottest atoms to escape. This cooling process can 
be continued until most of the thermally excited atoms are removed and, hence, the system can be considered as 
a condensate at T = 0. Differently from helium droplets, the properties of trapped gases can also be studied as 
a function of T, just stopping the evaporative cooling at the desired temperature. In this case, the gas has two 



components, a condensate and a thermal cloud, and the Gross-Pitacvskii theory of section III has to be generalized 



to properly include the thermal component. In the following, however, we restrict the analysis to T = |37| ]. 

B. Ground state 

Let us consider a typical helium cluster of 10^-10^ atoms. In its ground state, it is a self-bound liquid droplet with 
a rather flat density distribution, the central density being close to the one of the uniform liquid in the limit of zero 
pressure, i.e., n — 0.022 A~'^ (see Fig.l). Assuming the droplet to have constant density and sharp surface, its radius 
would be i? = tqN^^^ with rg ~ 2.2 A, which gives about 20 to 50 A, significantly larger than the average atomic 
distance. Actually, the true density profile has a rather smooth surface due to the large zero point motion of the 
atoms. The density decreases from the bulk value to zero within about twice the average atomic distance, i.e., about 
6 to 9 A depending on the precise definition of the surface thickness (see [Q and references therein) . The thickness 
of the surface region is comparable to the droplet radius for N of the order of 100 or less, but even for larger droplets 
the diffuseness of the surface can play a significant role. For instance, it makes the average density of the droplet 
lower than the bulk liquid value and hence the atoms are less bound to the droplet than to the uniform liquid. As 
already said, one can reasonably use the liquid drop model to express the energy per particle as a function of N: 

^=a, + asN-^/^ + acN-^/'' + ... (23) 

where ay ~ —7.2 K is the volume coefficient, i.e. the energy per particle in the bulk liquid (in the limit of zero 
temperature and zero pressure), and the next terms are surface and curvature energies. This type of formulae are 
frequently used to fit the numerical results of ground state calculations. 

The density distribution of a pure helium droplet is determined by the balance of the kinetic energy, associated with 
the zero-point motion of the atoms, and the potential energy due to atom- atom interaction. One can also produce 
density modulations in the droplets by picking up a foreign atom or molecule, acting as an external potential in which 
helium atoms readapt. This is, however, a local effect, in the sense that the main changes in the density distribution 
are limited to the first layers of atoms around the impurity. 

Differently from helium droplets, trapped condensates are not self-bound. They are produced and maintained in 
an external confining harmonic potential, of the form T4xt = (l/2)'7iu;^|^r'^ 139]. If this potential is switched-off, the 
condensate freely expands. 

The density distribution of a trapped gas is determined by the balance of three energy contributions given in 
Eq. (|lj). Since the sign of the interaction energy is fixed by the scattering length a, one may have either a repulsive 
or attractive mean-field potential, for positive and negative a, respectively. The case of negative a is very interesting 
for several reasons but it is not suitable for a comparison with helium. In fact, a uniform gas of atoms with negative a is 
unstable; it has negative pressure and it lowers its energy by collapsing into separate regions of high density. In a trap, 
quantum pressure can prevent this collapse, at least below some critical density p4[| , but such metastable condensates 
have no counterpart in the uniform limit, differently from the case of liquid helium and of trapped condensates with 
positive a. 

If the atoms were not interacting (a — 0) the condensate would have the form of the lowest single-particle state in 
the harmonic potential, i.e, a Gaussian having a width aho = [?i/(m'^ho)]^ • In actual condensates, with a > 0, the 
repulsion between atoms tends to lower the central density, thus increasing the condensate width. The equilibrium 
configuration is again fixed by the balance between i?qp, -Bext and £^int- It turns out that for Na/aho 3> 1 (Thomas- 
Fermi limit) the density in the central part is so smooth that the quantum pressure is negligible compared to both 
Eext and i^int- In this case, the GP equation ( pi] ) simply gives n{r) = g~^[^ — Vext('')], that is, a density having the 
form of an inverted parabola. In the same limit, one can easily calculate the energy per particle. 
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which contains both the harmonic potential parameters and the scattering length. Only in a narrow region at the 
surface, where the density vanishes, the interaction term becomes small and the density profile is determined by 
the balance between quantum pressure and external potential pO|; it deviates from the inverted parabola and goes 



smoothly to zero. This surface region gives a logarithmic correction to the Thomas- Fermi energy (p4|), of the form 
E/N ~ N~'^'^\aN . An example is shown in Fig.l, where one can see how the external potential affect the overall 
shape of the condensate (right panel), in contrast with the self-bound helium droplet (left panel). 

At present, several types of condensates are available in different labs (different geometries and atomic species). 
Let us consider a common situation of a condensate made by 10^-10^ atoms of sodium or rubidium, in a trap with 
cjho/(27r) ^ 10-100 Hz. The width of the condensate can be several /im, but wider condensates, up to 300//m, have 
been obtained in anisotropic traps. The scattering length a is of the order of a few nm and the ratio a/aho is about 
10"'^, so that Na/uho 3> 1; thus the density profile is very well represented by an inverted parabola with a thin 
surface region. The size of the system is much larger than the average interatomic distance, of the order of 10~^ m, 
as happens for large helium droplets. The central density is of the order of lO^'^-lO^^ cm~^, so that na^ is less than 
10"'^, which ensures the applicability of Gross-Pitacvskii theory. 

The relevant lengthscales are summarized in Table I for a ^He droplet with 10^ atoms and a trapped condensate of 
10^ atoms of *^Rb. One notes that the only difference in d, ^ and i? is a common rescaling, from nm (helium) to /im 
(trapped BEG). This is not true for a, which is of the same order of d in the dense fluid and much smaller than d in 
the dilute gas. 

A direct consequence of the very different values of na^ for the two systems in Fig.l is that the total density plotted 
in the right panel can be also identified with the condensate density, since quantum depletion is negligible (less than 
1%), while the same is not true for helium. In the latter case, the condensate density is expected to be about 10% of 
the total density in the inner part of the droplet, increasing to a maximum value at the surface where the system is 
more dilute mn (see also pi for a recent discussion). 

C. Excitations 

A peculiar feature of Bose superfluids is that their excitations at low energy correspond to collective modes, which 
can be described as fluctuations of the order parameter. This happens when their wavelength is larger than the 
healing length. 

For uniform and dilute gases, the spectrum of excitations is given by the Bogoliubov result (|l6|), which is phononic 
at low q and single-particle at high q. Low q excitations arc phonons even in liquid helium, with sound velocity c ~ 238 
m/s. However, since the fluid is dense and highly correlated, the interpolation between the phononic and single-particle 
regimes is more subtle. The phonon branch reaches a maximum at q ^ lA^^ and then forms a rather deep minimum 
at g ~ 1.9A~^, with energy ~ 8.7 K. The excitations near this minimum, whose dispersion is approximately parabolic, 
are called rotons. Their wavelength is of the order of the interatomic distances and slightly larger than the healing 
length, so that they still have collective character, but are related, roughly speaking, to a tendency to local order on 
the atomic scale. 

Helium droplets and trapped BEG have a finite radius R. As a consequence, the spectrum is discrete and the 
excited states have to be classified according to the number of radial nodes, n^, and the angular momentum quanta 
{I and m, if the system is spherical). The discretization is particularly important for the lowest energy excitations, 
whose wavelength is comparable with R, and correspond to oscillations of the whole system. 

An example of collective mode is the state with n^ = 1 and I = m = 0, which is a purely compressional oscillation 
(monopole, or breathing mode). By keeping I = m = and increasing rir, one finds density oscillations in the radial 
direction having wavelength smaller than R. If rir is large, so that the wavelength becomes much smaller than R 
but still larger than ^ (see Table I), these modes can be thought as stationary states of bulk excitations (phonons 
and/or rotons) which propagate radially and reflect at the surface. For large R, these discrctized states approach 
a continuum and the spectrum becomes closer and closer to the one of a uniform system, namely, the phononic 
Bogoliubov dispersion in dilute gas and the phonon-roton branch in liquid helium. 

Another interesting class of excitations is the one with n^ = and I > 2. These modes correspond to shape 
oscillations, or surface excitations, because the motion is mainly concentrated near the surface, the latter being 
displaced in and out periodically. The lowest one is the quadrupole deformation, 1 — 2. By increasing I one gets 
shape oscillations with more nodes in the angular direction, that is, surface waves with shorter wavelength. In the 
absence of external confinement as in helium droplets, and for large R, these states approach the collective waves of 
the liquid- vapour interface (ripplons), having dispersion uj'^ — q^a/{mn), where a is the surface energy and q is the 
wavevector parallel to the surface. In the case of trapped gases the dispersion is instead given by the law w^ = qF/m, 
where F = muuhoR is the trapping harmonic force calculated at r = _R |42|| . 

How can these collective states be excited and observed? In trapped condensate the answer is simple: one can 
modulate the external potential in order to force oscillations in the system in a classical way. With this method, one 
has already observed monopole, dipole, quadrupole modes, as well as higher multipolarities p3|-E5t|. Their frequency is 



measured with great accuracy and the agreement with theory is very good. The theoretical work basically corresponds 
to the numerical solution of the Bogoliubov-like equations (|l4|)-(|l5|). In the Thomas-Fermi limit, Na/aho ^ 1, one 
can also obtain analytic results by linearizing the hydrodynamic equations (|l8|)-(|l9|) pq ]. In spherical traps one 
obtains the dispersion relation cu{nr, £) = Who(2ri^ + 2nr£ + Sn^ + ^)^ ", which can be compared with the spectrum of 
noninteracting particles in the same harmonic potential, U!{nr,£) — uj^oi^rir + I) [ pJ[ . 

Phonon-like excitations with wavelength smaller than the condensate size can also be produced. For instance, one 
can suddenly switch-on a narrow laser beam, focused in the center of the trap. An optical dipole force acts on the 
atoms generating a wavepacket of excitations which then moves throughout the condensate as a sound wave. The 
velocity of this sound wave has been measured in E8|, finding good agreement with the prediction of Bogoliubov 
theory. Phonon-like excitations have been also generated in light-scattering experiments |4^ . 

Driving collective excitations in helium droplets, in a controllable way, is much less simple, since there is no confining 
potential to play with. However, one can produce and probe excitations by picking up impurities. A foreign atom or 
molecule, laying inside the droplet or at the surface, can be excited with a laser and this excitation may couple to the 
internal degrees of freedom of the droplet. This mechanism is at the basis of the so-called depletion method: a photon 
absorbed by an embedded impurity causes a transfer of energy from the impurity to the droplet via the creation of 
excitations. The latter eventually exchange their energy with atoms at the surface, letting them to evaporate through 
a quantum evaporation process. Thus, by monitoring the droplet beam with a mass spectrometer, one detects the 
spectroscopic transitions of the impurity which, in turn, contain information about the droplet itself pO| . This process 
has some interesting analogies with the quantum evaporation experiments in bulk helium ||51| , where the initial bunch 
of excitations (phonons and/or rotons) is produced by a heater immersed in the liquid, while the evaporated atoms are 
collected at a bolometer above the surface. There are also significant analogies with the sound propagation experiment 
in a trapped condensate B3|, where the "heater" is represented by the laser beam. In this case, if the trap is kept on, 
phonons do not evaporate atoms, but they may exchange energy with single-particle states in the outer part of the 
condensate. 

An example of how the superfluid environment affects the spectrum of embedded impurities is the observation of 
the phonon wing in the electronic excitation spectrum of a glyoxal (C2H2O2) molecule in helium droplets with about 
5500 atoms p2|. The observed features were found to be consistent with the existence of a phonon-roton spectrum of 
excitations inside the liquid droplet, with no low-lying single-particle excitations, as expected for a Bose superfiuid. 
Theoretical calculations have also shown that, for droplets of this size, the discrctized spectrum of collective excitations 
closely resembles the continuum of phonons and rotons of bulk liquid |2l|,^,^,Q . 

Not all impurities are bound inside helium droplets; some prefer to stay at the surface. In this case, they can be 
sensitive to surface modes. For large droplets, theoretical calculations show that the spectrum of these excitations, 
with / > 2, is close to the one of collective waves on a planar helium- vacuum interface (ripplons) |5q|. This also ensures 
that the liquid drop model (LDM) is accurate enough for many purposes. For example, using the liquid drop model 
for the excited states, together with a statistical method for the evaporation rate, the typical temperature of helium 
clusters in a supersonic beam was calculated in Ref. |g^. The predicted values, for both "^He and ''He droplets, were 
later found to agree with experiments [p4| , so providing an indirect insight into the spectrum of surface excitations. 

The fact that in a uniform superfluid system single-particle excitations are not allowed to exist below the low-energy 
collective modes is at the basis of the well known "Landau criterion". This says that, in order to create excitations 
in a superfluid by moving an impurity, one has to overcome some critical velocity, which is flxed by the spectrum of 
the excited states of the system. This follows from the conservation of energy and momentum. In a uniform system 
the critical velocity is Vc — inm[e{p)/p], where e{p) is the energy of an excitations carrying momentum p pq]. In 
liquid helium such a critical velocity has been the object of a longstanding investigation, involving the structure of the 
phonon-roton branch and the nucleation and dynamics of quantized vortex lines and rings. Recently an experiment 
was realized to test the applicability of the Landau criterion to helium droplets |57| ] . By studying low energy collisions 
of ''He and '^He atoms on ''He droplets, the authors found that ^He has a signiflcant probability to fly throughout 
the droplet without dissipation of momentum and energy, and they interpreted this as a consequence of the Landau 
criterion. The existence of analog critical velocities in dilute condensates is also currently investigated. In Ref. |5^] 
an external potential (a laser beam) is used to produce a hole in a trapped Bose-Einstein condensate and this hole 
is moved back and forth at variable velocity, playing the role of a massive impurity. Theoretical interpretations in 
terms of the Landau criterion have been already presented in Refs. p9| , |60[ . 

D. Rotations and moment of inertia 

When a system is described by an order parameter of the form (g) , a crucial consequence is the irrotationality of 
the superfluid flow, whose velocity fleld is given by v = {h/m)'VS, the phase S playing the role of velocity potential. 
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This is a fundamental feature characterizing the superfluid motion, which apphes to both hquid hehum and dilute 
condensates. In this connection, it is worth recalling that, despite the fact that in a strongly correlated superfluid 
|$P differs from the density n, the current associated with the superfluid velocity v, at T = 0, is still given by j = nv. 
The irrotationality of the superfluid flow has important manifestations like, for instance, the peculiar behavior of the 
moment of inertia and the occurrence of quantized vortices. 

The moment of inertia can be defined as the response, O = limo^o(i2)/^i to a rotating field of the form 
^rot = —^Lz, where Lz is the third component of the angular momentum operator. By rewriting the equations 
of hydrodynamics (|l8|)-(|l9|) in the frame rotating with angular velocity O, one gets the irrotational moment of inertia 

e^({^)'e.,. (25) 

where Ong is the classical rigid value. This reduction of the moment of inertia compared to the rigid value was 



observed in liquid helium in the 60's, so providing independent measurements of superfluid density |61|. In both 
helium droplets and trapped condensates, however, the direct measurement of Lz is not feasible and, consequently, 
the determination of Q requires more subtle probes. 

A suitable solution for trapped condensates is provided by the fact that, if the trap is deformed, the quadrupole 
and rotational degrees of freedom are coupled. This is well understood by taking a trap with frequencies ujx t^ ^y 
and considering the exact commutation relation [H, Lz] — im{uj'^ — u>y)Q, which explicitly points out the link between 
the angular momentum operator Lz and the quadrupole operator Q = ^^ Xiiji. Since the quadrupole variable can be 
easily excited and imaged in these systems, one can use it to get information on O. A successful method consists in 
preparing a deformed condensate, whose shape at equilibrium is an ellipsoid with different major axis along x and y, 
and suddenly rotating the magnetic trap by a certain angle around z. The condensate then responds with an oscillation 
of its major axis in the xy-p\ane around the new equilibrium orientation. This is called scissors mode. Using the 
hydrodynamic equations for an irrotational superfluid, one predicts the oscillation frequency uj ~ [cj^ + ^V^^"^ [ p2[ , 
which is different from the result of a classical gas, lo =\ ujxiiOy \. The explicit relationship between the scissors mode 
and the moment of inertia was derived in ||63| . Recent experiments |6j] have nicely confirmed these predictions. The 
observed moment of inertia clearly points out the superfluid nature of these condensates. 

In helium droplets, rotations can be induced and probed by acting on embedded impurities, namely, by studying 
their rotational transitions. Molecules like SFg, OCS, HCN, and many others have been already used (see Ref. [pSJ for 
a recent review). Their spectrum inside the helium droplet exhibits sharp rotational lines, as those of free molecules, 
except for a different value of their moment of inertia. The sharpness of the spectrum is already a signature that 
the cold matrix around the molecule is superfluid. A nonsuperfluid matrix would smear the spectrum so much that 
it could not be resolved; a demonstration of this effect has been given in Ref. [Q by replacing ''He (super) with a 
variable amount of ^He (normal). In this case, the '^He component provides the low- lying single-particle (particle-hole) 
excitations which cause the broadening of the rotational spectrum of the molecule. On the other hand, the fact that 
the observed moment of inertia is systematically larger than the one of the free molecule (up to a factor 4 or 5 for 
some of them) is an interesting consequence of a dynamical coupling between the rotational degrees of freedom of the 
molecule and of the superfluid environment. Ma ny exp erimental results are now available, but there are still open 
questions in their theoretical interpretation [p3| , p5[|67| , |6q| . One of the main problems is that the coupling between the 
rotor and the superfluid occurs on a microscopic lengthscale, comparable with both the interatomic distances and 
the healing length. On this scale, one is not fully ensured about the applicability of the superfluid hydrodynamic 
equations and/or of definitions of "local" superfiuid and normal densities. 

E. Rotations and quantized vortices 

Quantized vortices are striking manifestations of superfluidity. In liquid helium they have been the object of a 



longstanding investigation, starting from the pioneering ideas of Onsager and Feynman [ p9[ (sec |70| for a detailed 
overview). As soon as BEG was observed in trapped gases, the search for quantized vortices in these new systems 
became one of the primary goals. Recently, single vortices and arrays of vortices have been obtained with different 



techniques: by "phase imprinting" |71|, by using a laser beam "stirrer" |72(| , by rotating the magnetic trap ||73|| , and 
by rotating the thermal cloud during the evaporative cooling process WW- Conversely, quantized vortices have not 
yet been observed in helium nanodroplets. 

The quantization of the circulation in a superfluid directly follows from the assumption that the order parameter m) 
must be single- valued. Hence, the increment of the phase S over any closed path must be zero or an integral multiple 
of 27r. If the superfluid fills a simply-connected region, the irrotationality of the velocity field (V x v = 0) strictly 
implies that the increment is zero, and the superfluid can not rotate at all. In certain conditions, however, the fluid 
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may prefer to develop vortical lines, where the superfluid density vanishes. The increment of S around these lines can 
be an integral multiple of 27r and this implies that the circulation ^ v • dl is an integral multiple of the quantum h/m. 
The density distribution and the velocity field around a single rectilinear vortex line with circulation h/m can 
easily be calculated in the case of dilute gases by using the Gross-Pitaevskii theory. In fact, a quantized vortex can 
be regarded as a stationary solution of Eq. (noh of the form 
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where if is the azimuthal angle, r± is the distance from the z-axis, and ^y is a real function obeying the equation 
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The velocity field associated with the order parameter (^) takes the form v = {h/m)'Vip — z x r/r^ where z is 
the unit vector along z. It satisfies the irrotationality constraint everywhere except along the vortical line and gives 
rise to a total angular momentum given by Lz = Nh. The centrifugal term ex l/r'j_ in Eq. ( p7| ) originates from the 
peculiar behaviour of the velocity field of the vortical configuration. This term is responsible for the vanishing of the 
condensate density (j)^ along the z-axis. The "core" of the vortex, i.e., the region where the density is depleted, has 
a radius of the order of the healing length. The density profile and the energy of the vortex can be calculated by 
numerically solving the GP equation as done, for instance, by us in Rcf. U^. The OF theory can also be used to 
investigate the nucleation mechanisms and the stability of the vortex configurations (see | |7q ] for a recent review) ; the 
theory describes the main features of the vortices observed in trapped BEG and several predictions agree well with 
the available experimental results. 

What makes relatively easy to produce and manipulate vortices in trapped condensate is the possibility to induce 
rotations in a controllable way by means of laser beams or rotating traps. The vortical lines can be directly observed 
as "holes" in the density distribution. In fact, although the healing length is usually smaller than the resolution of 
the optical imaging devices, one can switch-off the confining potential, letting the condensate to expand, and wait 
till the expanding core is large enough to be observed [T^]. Alternatively, one can use the lowest collective modes of 
the condensate in the trap as a probe of vorticity. For instance, a vortical flow induces a frequency splitting of the 
two modes with / = 2 and m — ±2, which results in a slow precession of the quadrupole shape deformations [|78[ , 
which is analog to the splitting of the lowest circularly polarized modes of the vibrating wire in Vinen's "milestone" 
experiment with superfluid helium |79|. This method has been successfully used to estimate the angular momentum 
of vortical configuration [pOl, as well as to detect vortices which are not directly observable because they are bent 
or tilted |^l|. Finally, the 27r phase around a vortex can be also visualized in the form of a dislocation in the fringe 
pattern produced by the interference of two expanding condensates containing vortex lines 183] . 

In helium droplets there is no external field to play with and hence vortical configurations are not easy to produce 
and detect. Vortices might nucleate in the initial stages of the jet expansion, where droplets form and cool down, or 
in the collisions with molecules in the pick-up scattering chamber; but no quantitative predictions have been given 
so far about these processes. The stability of a vortex line in a droplet is also an open question. On the basis of 
purely energetic arguments, a vortex line is not stable since its formation costs a large amount of kinetic energy of 
the superfluid flow, and the ground state of the droplets is hence always vortex- free [p3| . On the one hand, however, 
if a vortex appears as a metastable state in some dynamical process, then one has to consider also the conservation 
of both energy and angular momentum in any possible decay mechanism, and one might discover that this state is 
indeed long-living, despite its large energy, as it has been shown for vortex configurations in trapped condensates. 
On the other hand, one can use dopant atoms or molecules to pin and stabilize the vortex, as the density functional 
calculations of Ref. p4] suggest. In this case, the dopant can be used also as a probe of vorticity, through the possible 
effects of the vortex flow on its rotational spectrum. An interesting perspective in this direction is the possibility to 
accommodate in the droplet a long chain of linear molecules, such as HCN, as done in Ref. [p5|. Such chains are 
expected to be suitable for the stabilization of a vortex line and, moreover, they can reveal in their spectrum the 
presence of the vortex flow, like for a vibrating wire in a nanoscale Vinen-like experiment. 



VI. CONCLUSIONS 

In this paper we reviewed some properties of both helium nanodroplets and dilute trapped condensates. Our purpose 
was not to give a systematic and detailed account of the available theories and experiments, but rather to discuss 
some relevant features which allow one to build a bridge between the two systems. We emphasized the existence of 
a common theoretical background and of many analogies in the measured properties, as well as of several interesting 
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differences. We mainly used the language of Gross-Pitaevskii and density functional theory because, on the one hand, 
it represents a very reliable mean-field approach for dilute gases and, on the other hand, density functional approaches 
has been successfully used also in helium droplets. Of course, the close connections between trapped dilute gases and 
helium droplets would equally emerge when examined from the viewpoint of more microscopic theories or stochastic 
Monte Carlo techniques. 

Once this bridge is established, one may include most of the recent advances in both fields within a unified per- 
spective, namely, the road towards a deeper understanding of the physics of interacting quantum fiuids. 
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N ^ 10^ 

t^ho = 27r(100Hz) 
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Table I. Approximate lengthscales (in units of m) for a typical helium nanodroplet and a spherical trapped BEC. 
a: range of interaction (helium) and scattering length (trapped BEC). d: average interatomic distance. ^: healing 
length, as in Eq. (p7[). R: radius of the system. The corresponding density profiles are shown in Fig. H. The energy 
per particle is about 6 K for the helium droplet and 60 nK for the condensate; the last value should be compared 
with huj\io ^ 5 nK. 
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FIG. 1. Left: density profile at T = of a He droplet obtained by solving Eq. ( |2l[ ) with the Orsay-Paris functional of 
Ref. [Eq]. Right: density profile at T = of a condensate of ^'^Rb in a spherical trap with Oho ~ 1-08 x 10-® m obtained by 
solving the GP equation ([ll[). 



15 



